Photon Recoil Momentum in Dispersive Media 
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A systematic shift of the photon recoil due to the index of refraction of a dilute gas of atoms has been observed. 
The recoil frequency was determined with a two-pulse light grating interferometer using near-resonant laser 
light. The results show that the recoil momentum of atoms caused by the absorption of a photon is nhk, where 
n is the index of refraction of the gas and k is the vacuum wavevector of the photon. This systematic effect must 
be accounted for in high-precision atom interferometry with light gratings. 
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The momentum of a photon in a dispersive medium is of 
conceptional and practical importance 13111111110]. When 
a photon enters a medium with index of refraction n, the elec- 
tromagnetic momentum changes from hk to Tik / n, where k is 
the vacuum wavevector of the photon. Momentum conserva- 
tion requires that the medium has now a mechanical momen- 
tum (1 — l/n)hk which for n > 1, is parallel to the propaga- 
tion of the light. The particles in the medium are accelerated 
by the leading edge of the light pulse and decelerated by the 
trailing edge. As a result, no motion is left in the medium 
after the pulse has passed. An absorbing surface without re- 
flection would therefore receive a momentum of hk per inci- 
dent photon. In contrast, a reflecting surface will recoil with 
a momentum of 2nhk per photon. In this case, the standing 
wave formed by the incident and reflected light pulse trans- 
fers momentum to the medium which remains even after the 
light pulse has left. This modification of the recoil momentum 
has so far only been observed for light being reflected from a 
mirror immersed in a liquid 00]. 

Recently, there have been discussions about what happens 
to an atom when it absorbs a photon within the medium. If 
one assumes that after absorbing the photon, no motion is left 
in the medium, then the recoil momentum should be hk fcj]. 
The same conclusion is reached when one assumes a very di- 
lute, dispersive medium with the absorbing atom localized in 
the vacuum space between the particles of the medium 11(11 . 
However, the correct answer is that the atom will recoil with 
a momentum of nhk, which requires particles in the medium 
to receive a backward momentum (for n > 1) due to the in- 
teraction of the oscillating dipole moments of the particles in 
the dispersive medium and the absorbing atom. So both for 
reflection and absorption by an atom, a photon in a dispersive 
medium behaves as if it has a momentum of nhk. 

In this paper, we study the recoil momentum of a photon in 
a dilute gas and show that its value is nhk. This has important 
consequences for atom interferometers using optical waves to 
manipulate atoms by the transfer of recoil momentum. High 
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precision measurements of the photon recoil are used to de- 
termine the fine structure constant a 

Further improvements in the accuracy of photon recoil mea- 
surements, combined with the value of a derived from the 
(g-2) measurements for the electron and positron B17U18U19I1 . 
would provide a fundamental test of QED. The accuracy of 
the best photon recoil measurements are limited by the uncer- 
tainty in the correction to the photon recoil due to the index of 
refraction. At low atomic densities, where atom interferome- 
ters usually operate, the index of refraction effect is relatively 
small. Here we operate an atom interferometer with Bose- 
Einstein condensates, which have much higher density than 
laser cooled atomic clouds, and observe how the index of re- 
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fraction modifies the atomic recoil frequency hu rec = 
where p and m are the atomic recoil momentum and mass re- 
spectively. 

The essential idea of our experiment is to measure the re- 
coil frequency interferometrically using a two-pulse Ramsey 
interferometer. The two pulses are optical standing waves sep- 
arated by a delay time r (Fig. [I]). The first pulse diffracts the 
atoms in an 87 Rb condensate into discrete momentum states. 
During the delay time r the phase of each momentum state 
evolves at a different rate according to its recoil energy. The 
second pulse recombines the atoms with the initial conden- 
sate. The recombined components have differing phases lead- 
ing to interference fringes that oscillate at the two-photon re- 
coil frequency. By measuring the resulting frequency as a 
function of the standing wave detuning from the atomic reso- 
nance, we found a distinctive dispersive shape for u rec that fit 
the recoil momentum as nhk, where n is the index of refrac- 
tion of the gas. 

The experiment was performed using an elongated 87 Rb 
Bose-Einstein condensate created in a cloverleaf-type Ioffe- 
Pritchard magnetic trap using the procedure described in 
Ref.|20]. The condensate, containing 1.5 x 10 6 atoms, was 
produced in the |5 2 5 , i/ 2 , F = 1, n%F = —1) state, and had a 
Thomas-Fermi radius of 8 (90) fim in the radial (axial) direc- 
tion, and the magnetic trap had a radial (axial) trap frequency 
of 81 (7) Hz. 

The BEC was illuminated from the side with an optical 
standing wave created by a retro-reflected, 7r-polarized laser 



2 




(a) 



2nlik 



Onfik 



-2nfik 



Release from Trap 
Image inTOF 



FIG. 1 : Kapitza-Dirac interferometer. The first pulse outcoupled a 
small fraction of atoms into the | ± 2nTik) momentum states. The 
outcoupled atoms moved within the initial condensate. After a vari- 
able delay r a second pulse was applied and outcoupled atoms that 
interfered with those outcoupled by the first pulse. The laser beam 
was applied perpendicular to the long axis of the condensate; the po- 
larization, E, was parallel to it and to the applied magnetic field bias, 
B. 600 ijls after the first pulse was applied the atoms were released 
from the magnetic trap and imaged after 38 ms of ballistic expansion. 
The field of view is 0.5 mm x 1.5 mm 



beam. The polarization of the beam was optimized by sup- 
pressing Rayleigh superradiance 1 20] . The laser was detuned 
from the 5 2 S 1/2 ,F = 1 -> 5 2 P 3/2 ,F = 1 transition at 
A = 780 nm, and had a linewidth 7 much smaller than T, the 
natural linewidth of the transition. The intensity of the 5 /is 
long pulse was set to outcouple « 5% of the atoms i nto each 
of the \±2nhk) momentum states with no appreciable popula- 
tion in higher momentum states. This ensured that the density 
of the original condensate was nearly constant throughout the 
measurement. After a variable time r, a second identical pulse 
was applied. The time between the first pulse and the shutoff 
of the magnetic trap was fixed at 600 /is, which was less than a 
quarter of the radial trap period. The momentum distribution 
of the condensate was imaged after 38 ms of ballistic expan- 
sion, long enough for the momentum states to be resolved. 
The images were obtained using resonant absorption imaging 
after first optically pumping the atoms to the 5 2 £1/2, F = 2 
state. To compensate for spontaneous light scattering from 
the standing wave, the density of the condensate (and associ- 
ated mean field shift) was determined by applying a single 5 
lis pulse to the condensate, and then immediately releasing it 
from the magnetic trap. The number of atoms in the conden- 
sate was determined by integrating the optical density of the 
absorption image. The optical density was calibrated by fit- 
ting the Thomas-Fermi radius of unperturbed condensates in 
time-of-flight jH. 

The recoil frequency was found by fitting the oscillations 
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FIG. 2: Interference fringes oscillating at the recoil frequency, (a) 
Absorption images for r = 10 — 50 /jls. The detuning was A/2tv = 
+520 MHz relative to the F=l — > F'=l transition. The field of view 
is 0.5 mm x 1.5 mm (b) Fraction of atoms in the \0hk) momentum 
state as a function of r. The fringes were fit using Eq. Q. The fitted 
frequency was uj — 2tt x 156268(39) Hz and the decay constant 
t c — 461(25) /is. The signal was normalized using the total atom 
number in all momentum states. 



in the fraction of atoms in the \0nhk) momentum state for r 
ranging from 10 — 590 fis (Fig.|3. The interference fringes 
were fit to a cosine function with a gaussian envelope: 



Aexp ( 1 cos(4cjt 



C. 



(1) 



The observation of up to ten oscillations provided a precise 
value of the recoil frequency. The origin of the damping time 
r c and of the offset C will be discussed later. 

Figure[3] shows our measured values for uj/2tt as a function 
of the detuning, A/27T, relative to the F=l — > F'=l transition. 
The measured values for the recoil frequency clearly follow 
the dispersive shape of the index of refraction, and they are 
in good agreement with the hypothesis that the recoil momen- 
tum is proportional to the index of refraction. The variation in 
uj/2tt as a function of the detuning was 2 kHz across the reso- 
nance, much larger than the statistical error on the frequency 
fits of less than 100 Hz. This conclusively shows that the mo- 
mentum transferred to the atom when a photon is absorbed is 
nhk. 

We now discuss in more detail how the atoms interact with 
optical standing waves. For the short duration of the applied 
pulses (5 /is) we can assume that the atoms do not move 
during the pulse and ignore the kinetic energy of the atoms 
(Raman-Nath approximation). The interaction can then be de- 
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FIG. 3: Recoil frequency as a function of detuning A/2tv showing 
the dispersive effect of the index of refraction. The average density of 
the condensate for the solid points was 1.14(4) x 10 14 cm -3 , giving 
rise to a mean-field shift of 880 Hz. The shaded area gives the ex- 
pected recoil frequency including the uncertainty in the density. The 
dashed line is at u — An 2 uj re c + pU/h, the expected value without 
index of refraction effects. The data shown as open diamonds had 
increased spontaneous light scattering due to cr^ light contamination 
in the laser beam. The increased light scattering led to a lower initial 
density in the condensate, thus leading to a smaller mean-field shift 
and index of refraction. The contamination allowed Awif = ±1 
transitions, thus for small detunings the proximity to the |1, — 1) — > 
|0 / , 0) transition located at A/2tt = —72 MHz resulted in higher 
spontaneous scattering rates. The open points have been scaled up- 
ward to correct for this lower density. 



scribed by the application of the AC Stark potential due to the 
standing wave 



V(z) = ^^sm 2 (nkz), 



|^(r)> = \ijj ) ± 2nhk)e- i4n2uj ^ T + J o (0)|Onftfc>) . 

(3) 

At t = r a second pulse is applied that partially recombines 
the momentum states. After applying the two pulses, the prob- 
ability of finding the atoms in the the \0nhk) momentum state, 
Po = KVK 7 " + t p )\0nhk)\ 2 is given by 

po = 4(0) + 4 [J 2 (0) 4(0) + 4(0)] cos(4nV ec r). (4) 

As a function of r the density of the zero momentum peak 
oscillates at An 2 uo rec . 

So far, we have ignored the motion of the atoms during 
the delay time r. However the amplitudes of the recombined 
components will only interfere where they spatially overlap. 
After the first pulse, the atoms in the | ± 2nhk) momentum 
states move within the initial condensate with the recoil ve- 
locity (~ 12 /im/ms). As the overlap between the recoiling 
atoms and those at rest decreases, the interference fringes de- 
cay. The overlap integral for this decay can be approximated 
as a Gaussian with time constant, r c « 0.75 Rtf/v^c* where 
Rt f is the Thomas-Fermi radius of the condensate, and v rec 
is the recoil velocity [24]. 

The index of refraction for the condensate is derived from 
its macroscopic polarization P. For a two level system, 

,2 



P = x^oE = ip 
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where x i s me atomic susceptibility, e Q is the permittiv- 
ity of free space, p is the dipole matrix element, p is 
the atomic density of the condensate. In this experiment 
the light was 7r-polarized and detuned by A from the 
5 2 Si /2 ,F = 1 -> 5 2 P 3/2 ,F / = 1 transition. For this 
polarization the selection rule is Amp = 0, and there are 
two allowed transitions from \F = 1, = —1) — > — 1) 
and 1 1,-1) — > |2',-1) which are separated by 157 MHz. 
Including both transitions in the derivation, the index of 
refraction, n = y/1 + x, is given by: 



where A is the detuning between the optical frequency and the 
atomic transition, and ujr is the Rabi frequency. This equation 
is valid for large detuning, A 2 ^> T 2 /4, where Y is the natural 
linewidth of the transition. The short pulse limit, describing 
Kapitza-Dirac scattering, is valid for short interaction times 
t p , where t p <C l/u; rec w 40/iS. The first pulse outcouples a 
fraction of atoms into the momentum states | ± 2£nhk), where 
the population in the I th momentum state is given by Pi = 

4(0) ESQ, where for a square pulse, = and J £ 

is the £^-order Bessel Function of the first kind. For < 
1 a negligible fraction of atoms is diffracted into states with 
t > 1, and we can restrict our discussion to the | ± 2nfik) 
momentum states. For our experimental parameters 6 — 0.45. 
During the delay time r the phase of the | ± 2nhk) states 
evolves at a faster rate than the \0nhk) state due to the recoil 
energy, E rec = 4n 2 h(j rec , hence the wavefunction will evolve 
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where Ai and A 2 are the detunings relative to the F = 1 — > 
F' = 1 and F 1 — 2 transitions, respectively. This equation is 
valid in the limit pX 3 <C 1 [25], where A = X/2tt. For our 
experimental parameters pX 3 w 0.2. In addition to the index 
of refraction shift, the observed recoil frequency has a mean 
field shift [26]; the atoms in the | ± 2nhk) state have twice the 
mean field energy of those at rest due to the exchange term 
in the interatomic potential. Including both the mean field 
shift and the index of refraction, the frequency of the observed 
interference fringes should be: 



uj = 4n uj r 



pU 

h 



(7) 
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where pU = 47rh 2 ap/m, and a is the s-wave scattering 
length. The density p = 4/7p G , where p Q is the peak conden- 
sate density and the factor of 4/7 is due to the inhomogeneous 
condensate density. 

When the interference fringes were fit using Eq. CJ, the 
average values for the amplitude A and offset C for all of the 
data points were 0.12(3) and 0.82(4), respectively. This is in 
reasonable agreement with the expected values of A = 0.18 
and C = 0.81 for = 0.45. For a Thomas-Fermi radius of 
8 pm we would expect a decay time r c « 500 ps. There 
was an unexplained shift in the fitted value for the decay time 
r c between the red and blue side of the resonances, on the 
red side the average value was 347(20) ps and on the blue, 
455(40) ps. 

The quadratic dependence in Eq. on the index of refrac- 
tion can be understood by considering the diffraction of atoms 
from the light grating created by the standing wave. When the 
first pulse is applied, the standing wave creates a grating with 
a periodicity d = A/2. Within the condensate the index of re- 
fraction modifies the grating period by n, since A' = A/n. The 
recoil momentum of atoms that diffract off the grating will be 
changed by ±2hk f , again within the condensate k' = nk and 
the velocity of the atoms is modified. Assuming n > 1, when 
the second grating is applied the atoms have moved farther by 
a factor of n and the grating is shorter by n, changing the time 
scale for the interference fringes by a factor of n 2 . 

The increase in the momentum transferred to the atoms can 
also be explained by considering the momentum transferred to 
atoms by a classical field. A derivation using the Lorentz force 
applied to the atoms during the absorption of a photon can be 
found in [6]. In a dielectric medium with n > 1, the magnetic 
field and therefore the Lorentz force are not modified. How- 
ever, the electric field is weaker and therefore it takes longer 
for the atom to perform half a Rabi cycle and be transferred 
to the excited state. During that longer time, the Lorentz force 
imparts a momentum to the atom which is larger than ftk. 

For Kapitza-Dirac scattering, atoms are diffracted symmet- 
rically into the into | ± 2tnhk) momentum states, so mo- 
mentum is clearly conserved. However for processes such as 
Bragg scattering, where the atoms are scattered in only one 
direction the index of refraction has an additional effect. As- 
suming a 7r/2 pulse with counter-propagating beams, where 
half the atoms are diffracted, for n > 1 the recoil momen- 
tum is a factor of n higher than the vacuum momentum. For 
momentum to be conserved, the remaining atoms must recoil 
backwards with momentum p = 2(n — l)£hk. For small frac- 
tional outcoupling the effect is negligible, since the extra mo- 
mentum is distributed among the remaining condensate. How- 
ever if a large fraction of the condensate is outcoupled and t 
is large, this effect could potentially be resolved in ballistic 
expansion. 

We have discussed here the dispersive effect on the photon 
momentum near a one-photon resonance. An analogous effect 
occurs near two-photon resonances. In this case, the atomic 
polarizability is determined in third-order perturbation theory, 
and the resulting index of refraction has a sharp, narrow dis- 



persive feature near the two-photon resonance [27]. In recent 
experiments at Stanford 111 311 . such two-photon effects have 
been the leading source of uncertainty in high-precision deter- 
minations of atomic recoil frequencies and the fine- structure 
constant a. 

In conclusion we have measured a systematic shift in the 
photon recoil frequency due to the index of refraction of the 
condensate. This is the first direct observation of the atomic 
recoil momentum in dispersive media. For high atomic den- 
sities, this shift can have a significant effect on atom inter- 
ferometers, and is of particular importance for precision mea- 
surements oih/m and a with cold atoms H 1 3ll 1 6fl . 
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